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A generic outcome of theories with scalar-tensor coupling is the existence of infla- 
tionary attractors, either power-law or de Sitter. The fluctuations arising during 
this phase are Gaussian and their spectrum depends on the wavenumber k ac- 
cording to the power-law k l ^ l ~ p \ where p is the inflationary power-law exponent. 
We investigate to which extent these properties depend on the coupling function 
and on the potential. We find the class of models in which viable attractors 
exist. Within this class, we find that the cosmic expansion and the scaling of 
the fluctuation spectrum are independent of the coupling function. Further, the 
analytical solution of the Fokker-Planck equation shows that the deviations from 
Gaussianity are negligible. 
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1 Introduction 



Non-minimally coupled theories of gravity are models in which a direct coupling of the form 

1 



A NM c = J d 4 x ^g 



2 U(<P)R 
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between a matter field and the curvature scalar R is added to the gravity Action. The 
interesting history of such theories in general relativity and cosmology dates back to the 
Dirac suggestion and to Jordan-Brans-Dicke (JBD) models^. Terms with non-minimal cou- 
pling (NMC) arise quantizing fields in curved space!, in multi-dimensional theories! like 
superstrings or supersimmetry and in induced gravity theories!! On the other hand, La and 
Steinhardt's extended inflation^ (EI) and many variants within the context of 'first-order 
inflation' □ reconcile a first-order phase transition with inflation modulating the bubble nu- 
cleation rate. The cosmological kinematics of NMC theories is also discussed in Ref. (|7|). 

No known fundamental principle predicts the functional form of /(</>). One can find in 
literature many ansatze: f(<fr) = <ft 2 , as in the induced gravity model and some version of 
EI; f((f)) = a + bcj) 2 + c0 4 + as in the hyperextendedB design of EI; /(</>) = e D ^, as in 
dimensionally reduced Kaluza-Klein theories, and so on. In Ref. || , non-minimal derivative 
couplings in /(</>) have also been taken into account. Models with a generalized kinetic term 
like o;(0)0 ;A1 ;M /0 and a coupling term $i? have also been considered! (the case uj = const 
is indeed the JBD theory and the original formulation of EI); up to a redefinition of <p these 
models are equivalent to introducing a coupling like (|1]). Aim of this paper is to investigate in 
a systematic fashion the dynamical features of models with a general coupling /(</>), without 
confining ourselves to a specific choice, except for the following requests: that in the large 
limit the function /(</>) grows faster than 2 (the case /(0) ~ 2 being already well-known), 
and that the potential V(<fi), at least for large 0, can be written as 

V(d>) = A/ M (0) , (2) 

where A and M are arbitrary non-negative constants. [] We also assume V(<p), f(<f)) > 0, 
the equality being verified when eventually falls into its stable minimum. It is remarkable 
that, with the above assumptions, it is possible to find analytically, without specifying nor 
the coupling nor the potential, the class of inflationary attractors, the order-of-magnitude 
amplitude of the primordial perturbations and the solution of the Fokker-Planck equation for 
the stochastic fluctuations. Since /(0) is left largely unspecified, we speak here of theories 
with general coupling , instead of non-minimal coupling models. We also give examples of 
numerical phase-space (PS) of the investigated models. 

The plausibility of inflation depends on how large is the set of initial conditions which 
evolve to inflationary expansion of the scale factor: the existence of an inflationary attractor 
extending to values of the field as large as possible, up to the Planck boundary, is then a 
very desirable property of a cosmological model. A first result of the present paper is that 
we determine the class of models with coupling different from the standard choice <j) 2 R which 
have asymptotic inflationary attractors. We find that, in most models, either the Universe 
does not inflate at all, or it does so toward the wrong phase-space direction. The main 

a After this work was completed we learned of independent work of De Ritis and coworkers^, who 
obtained exact solutions in particular cases of the above functional relation. 
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result is that, in those cases in which a successful inflation is recovered, the kinematics of 
the Universe and the main features of the fluctuation spectrum do not depend on /(</>). 



2 Inflationary attractors and the perturbation spectrum 

We start from the following Action [we assume Planck units, G = c = % = 1 and signature 
(+ )1 



A = d A x y/=g 



107T L L 



(3) 



Great simplification is obtained passing to the the so-called Einstein frame, i.e. deriving the 
field equation in the new metric = e^g^ with 

2uj= ln|l- 7 /|, (4) 

where 7 = 8tt^ < 0. We will work here assuming £ < 0, which guarantees that the new 
metric is non-singular. The old metric will be referred to as the Jordan frame. The Einstein 
equations in the new frame are then 



G " U 1 - 7/ 



37 2 r 1 



(5) 



where is the usual energy-momentum tensor for the scalar field <fi. From now on, we 
confine ourselves to a homogeneous, isotropic and spatially flat metric with cosmic factor 
a(t) and Hubble function H(t) = a/a. Unless otherwise stated, all quantities are meant to 
be expressed in the rescaled variables: a dot denotes differentiation with respect to the new 
time t = J ^dt', the Hubble function H is in terms of the new cosmic factor, and so on. If 
Kfa, = (p-,^(p;u — \g l iv4'\a ( P' ,a denotes the kinetic sector of T^ u in the original frame, the kinetic 
sector in the new Einstein frame can be written in the form = F 2 (<f))K llu , where 

F 2 (0) = [16tt(1 - 7/) + 3 7 2 /' 2 ]/16vr(l - iff . (6) 

A canonical kinetic sector is then obtained defining a new field 

i>=\ . (7) 



This integral is not easily done, even with simple choices of the coupling f(<f>). However, 
a powerful simplification is attained in the limit \<p\ — > 00, the same limit in which the 
effect of the coupling is greater and inflationary attractors are found. Let us denote in the 
following with a subscript a quantity evaluated at the initial time. Then, if f' 2 grows 
faster than /, we can integrate Eq. (|7|) into ap = ln[/(0)//j] , and thus /(</>) = ftexpcip, 
where c = (167T/3) 1 / 2 . It then follows that the initial value of ip is ipi — 0. When f' 2 grows 
exactly like /, i.e. for / = (ft 2 , one has instead c(£) = 2[y/(6£ — 1)] 1//2 . 

Suppose now that the following relation holds between the coupling / and the potential: 
V = Xf M . This relation is verified, for instance, if both V and / are power functions of the 
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field, or if they are both exponentials. The potential in the Einstein frame, in the same limit 
as above, is then 

W) = A 7 - 2 r = /3exp(c/^), (8) 

where j3 = ft^l~ 2 > and p, = M — 2 . We then explore the bidimensional parameter space 
(£, M). Remarkably, an asymptotic flat potential able to drive a slow-rolling inflation in a 
FRW metric with scale factor a(t) is found in a single case: V ~ f 2 . In all other cases, we 
find in the if>,if> phase space the attractor trajectories 

^ = Ae c ^ /2 , A 2 = 2/V/(9 - p 2 ) , (9) 

where the negative root for A is to be chosen, corresponding to <ft decreasing toward its stable 
minimum, where the Friedmann behavior takes place. This corresponds to trajectories 

u = cA/r /2 & +i)/2 , (io) 

in the (4>, 0) plane. In the Appendix, solutions of this kind will be proven to be attractors; 
graphic evidence of their attractive properties is provided by our Poincare projected phase 
spaces (Figs. 1-4). A power-law expansion (for M^2) takes place on the attractors: 

a = a i (l + t/r) p , p = 3/p 2 (11) 

(in the conformally rescaled time), where r = —2/Acp > 0. The cancellation of the coupling 



constant £ from ( 11) occurs only in this class of models: in the ordinary coupling / = <p 2 
the cancellation does not take place. Notice that the solution (|ll|) as well as the attractor 
(|) above are exact solutions as long as we take the potential (§). From (|HD we see that 
H = Hi(l + t/r)- 1 , where Hi = p/r. 

The behavior of the cosmic scale factor is accelerated, hence inflationary, only if M < 
y/3 + 2 ~ 3.7. However, for M < 2 the model enters an eternally inflating phase, since then 
dU /dip < and the field if) is pushed to ever growing values, never reaching the Friedmann 
phase located at if) — > — oo [from the mapping if> — > (f> follows that if) — ► — oo when /(</>) 
vanishes]. The only allowed range is then 

2<M<3.7. (12) 

Outside this narrow range, the model studied here does not allow a successful chaotic infla- 
tion. When f = <f> 2 the expansion is modified to p = p(l — l/2£). This is inflationary for 
any M if £ < is small enough!. 

In the original Jordan frame the cosmic scale expansion is still a power-law in the Jordan 
time, but with an exponent pj = (p + /x)/(l + n), which is inflationary if p > 1. 

Since we have our theory in the Einsteinian form, we are allowed to employ the standard 
formalismliHl to derive the perturbation spectrum at the horizon-crossing epoch. We will see 
that we can determine the spectrum without specifying the coupling function. We assume 
here that the horizon-crossing condition is to be found in the Einstein frame, i.e. that the 
condition only involves the rescaled variables [see Ref. (|TT|)]. 

The rms amplitude 5p/p of the perturbations on a comoving scale k~ x at horizon crossing 
(after the end of inflation) is given by 

*P/P=(*V*) M , (13) 
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times an order-of-unity depending on the power-law exponent p. Along the solu- 

tion (0) we have 

H 2 /^ = {T/2)c^H 2 i {l + t/r)- 1 (14) 

Now, since k = aH = ajifj(l + t/r) p_1 , eliminating the time between the latter equation 
and fln| ) we obtain 

5p/p = (p/2)cfiH l (k/a t H i ) 1 ^ , (15) 

where the last factor expresses the scale dependence of the spectrum. Since p > 1, the 
spectrum grows with the scale, as it is commonly found in power-law inflation, and becomes 
flat for p — > oo (slow-rolling inflation). Our result shows that this conclusion does not depend 
on the choice of /(</>), as long as V = Xf M . 

3 Stochastic approach to fluctuations 

To compare the theory with the observed large-scale structure we need also the higher-order 
moments of the fluctuation distribution. In particular, it is astronomically interesting to 
detect possible deviations from Gaussianity of the primordial spectrum^. We investigate 
then the stochastic properties of adiabatic fluctuations during a power-law inflation driven 
by the JBD field 0. The pioneering work in this context is Ref. (O). Since we have 
now an Einsteinian theory with an exponential potentiate, we may analytically solve the 
Fokker-Planck equation (FPE) for the probability distribution of the fluctuations. The only 
other known cases of analytically solvable FPE are the minimally coupled [f{4>) = 0] case 
with quartic or exponential potentials^!. In power-law inflation, contrary to slow-rolling 
inflation, the term ip is not generally negligible in the field equation for ip. However, we 
show that its contribution is not dominant in the cases investigated here, at least for M 
not close to its upper bound 3.7. The condition \ip\ <C \U'\ gives in fact, when applied 
to the attractor solution (|j), A 2 /2f3 <C 1, that is, (c/z) 2 <C 24tt. We have, for instance, 
A 2 1 '2/3 ~ 0.12 for M = 3, which is moderately satisfactory, and A 2 /2(3 ~ 0.03 for M = 2.5, a 
good approximation. We can then safely put ip = —U'/3H and write the Langevin equation 
in the useful time variable a = J Hdr 

d T = ~^H + J^9Hv(a) = -f- + g^e^rjia), (16) 
da 3H Z V o7T 87r v 

where i](a) is the stochastic Gaussian noise due to the super- horizon stretching of high 
frequency modes, and g is an order-of-unity constant!^ which depends on the power-law 
exponent p. We have made use in flT6"|) of the (0, 0) Einstein equation in the Einsteinian 
frame in slow-rolling regime, 

3H 2 = 8vrA7" 2 /f exp(c/i^) . (17) 

We will denote here with a subscript 'cZ ' the unperturbed quantities, i.e. the solution to the 
zero-noise classical equations of motion, to distinguish them from the stochastic variables 
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which appear in the Langevin and Fokker-Planck equations. Let us define the new field 
variable 



dip 



27 



_J~^ e -c^/2 



from which Eq. (|T6|) for / — > 00 is simplified to 

d^/da = ty/p + 7)(a) . 
The FPE associated to Eq. (|19j) is in the Ornstein- UhlenbeclH form 
dW 



da 



1 <9( W) d 2 W 



d{D x W) + &(D X xW) 



dV 2 



(18) 



(19) 



(20) 



where Di = ty/p is the drift coefficient and Du = 1 is the diffusion coefficient. One sees that 
the definition domain of \I> is (0, +00) for /(</>) non-negative definite (remember that 7 < 0). 
Then we must impose the reflecting boundary conditions S(^f = 0) = , S(^/ — > 00) = 
where 



S = -D X W + d{D n W)/&& , 

is the probability flux. This introduces a 'reflecting' term in the solution^, which t 
to be negligible at large times, when < ^ > moves to infinity. The general solution^ 
with initial condition 0) = 5(^> — is then found to be 



(21) 

irns out 
lof (EDI) 



2np (e 2a /p - I) 



cxp 



'2p(e 2a /p - 1) 



+ exp 



2p(e 2a /p - V 



(22) 



where \l/ c / = ^ie a ' p . It is easy to see that all the moments of ip are defined. In fact, the 
integral 



mm =< r >-- 



ip K W{ip,a)dip. 



(23) 



where W(ip,a) = W(^/,a)\J\, and \J\ = d^/dip, converges for any k. The same is not 
always true for the moments of the original field 0. For instance, if f(4>) = <p 2m , we have 
that all the moments < <p K > for which k > mfi are divergent for <p — > 00. 

The distribution of the field ip is non-Gaussian. We are interested in estimating its 
deviation from Gaussianity through the skewness coefficient 



M 3 ^-ii cl )/Ml /2 (ij-ij cl 



(24) 



In the following we will neglect the contribution from the reflecting term, so that we have 
Mi(ip) = ip c i. We may first calculate approximately the second central moment of ip, i.e. the 
variance, in the limit of small \ip — ip c i\- We will see that this is completely justified in viable 
theories of inflation, due to the smallness of H, which sets the amplitude of the diffusion 
around the classical trajectory. Expanding — ty c i\ in powers of \ip — ip c i\ we can write the 
solution W(ip) to the first order as a Gaussian with mean ip d and variance 



D 



V9\e 2alv ~ 



l)[A 7 - 2 /fe c ^] 



(25) 
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We can now recognize in the square brackets, apart order-of-unity factors, the classical value 
of H 2 . It follows that, for large a, 

yf< ty-^ cl ) 2 > « H cl e a / p = H cl (t/r) . (26) 

For i> rwe have Hdt/r ~ Hi. In order to avoid exceedingly large perturbations on the 
microwave sky, Hi must be smaller than, roughly, 10~ 4 in Planck units, as it is evident 
from Eq. fllSP . It turns out that the deviation of if) from ip d is indeed very small, which 
justifies the approximations made above. We may now calculate the third central moment 
M 3 (if) — ipd) in the following way. Let us consider a Gaussian distribution VT(^) normalized 
to unity, with mean ty c i and variance a. We transform it under the mapping \I/ = ^(if>) to 
[\I/(-?/>)] | |, where ty' = d^/dip. The inverse mapping is if) = if>fif). The mean of if) is now 
V'rf = ipi^ci), while the variance and the higher-order moments have to be calculated. The 
new distribution is W(ip) = W[^/(ip)}\^/'\. In the limit \ip — ipci\ -> we may expand ^{if>) 
and ty' in Taylor series around ipd. It follows 

Wty)\V\ = (2nD)- 1 / 2 exp{-^ [ty - ^ cl f - Rty - ^)] } , (27) 

where D = cr|*' d |- 2 and R = 2D^/^' d . Here * d mean and respectively, 
calculated in if) — if> c i. The distribution ( p7|) is again a Gaussian curve, but with mean, 
variance and normalization different from the original ones. After some simple manipulations 
it is easy to find 



Maty -,!>*) = - [3D + — ) . (28) 




Let us specialize to our model. We have then that D is given in (P5|), while R = —cfiD. In 
the limit of small D, i.e. small H 2 , we have 

s = -3cfiD 1/2 /2 H d (t/r) . (29) 

The deviation from a Gaussian distribution turns out to be indeed negligible as long as the 
initial value for H is less than 10 -4 in Planck units. In other words, the need for a smooth 
microwave background has the consequence of a very small deviation from Gaussianity. 
Perturbations which crossed out the horizon when H ~ 1 can have a very skewed distribution, 
but they are now much larger than the observable Universe. Similar conclusions hold for the 
higher moments of the distribution. 

4 Conclusions 

One of the most appealing aspects of chaotic inflation is that even the simplest models possess 
the required properties for a successful inflation. Almost all the possible initial conditions 
for a scalar field dominated Universe with a simple potential lead to a conspicuos period of 
inflation. This feature implies, in the dynamical phase-space of such theories, the existence 
of attractor trajectories over which the other trajectories converge. It turns out, contrary to 
the uncoupled case, / = 0, and to the usual f = 4> 2 case, that models with generic coupling 



7 



satisfy the chaotic inflation requirements only when 2 < M < 3.7. The resulting power- 
law expansion is independent of /(</>). This power-law case is noticeable because it can be 
adopted as a chaotic hyperextended inflation with a kinematics independent of the coupling 
function. 

The scale-dependence, k l K l ~ p \ also independent of /(</>), is an interesting signature of 
power-law inflation, possibly testable with the data from COBE. Solving the Fokker-Planck 
equation we can evaluate higher-order moments of the fluctuation distribution. This shows 
that the deviation from Gaussianity on observable scales is negligible in all acceptable infla- 
tionary models of the kind studied here. 

Finally, we prove analytically in the Appendix the stability of the attractors. 



Appendix 

Numerical phase-spaces show with full evidence the stability of the trajectories labelled 
as A-type in the text. Here we prove it analytically for the case M > 2. As it was shown 
in Sect. 2, the phase-space in the large limit is governed, in the rescaled variables, by the 
dynamical system (x = if) , y = ip) 

x = y, y = -3H(x, y)y - U'{x) , (30) 

where U' = /3cfie Cflx . There is an attractor trajectory 

yo = Ae c ^'\ (31) 

with c/x > and A < 0. Here we assume for simplicity that H is calculated along the 
attractor H = [87r(l/2 + (3/A 2 )/3] 1 / 2 \yo\, freezing the x, ^-dependence of H on the assumed 
solution fl31~|) . Dealing with the full function H = H(x, y) only entangles the algebraic details 
without changing the conclusions. To show that ([H]) is indeed an attractor, let us define a 
new set of coordinates u, v, such that u = defines the trajectory (j3lD : u = y—Ae c ^ x ^ 2 , v = 
y + Ae^/ 2 . Under this mapping, the dynamical system ( |30"D becomes 

2ii = ^(u 2 -v 2 )-3H(u + v)-2U' , 

2v = ^V-u 2 ) -3H(u + v) -2U'. (32) 

We now linearize (|32]) around xq, yo = Ae cflXo ' 2 , i.e. around a generic point on Ae c ^ 2 . This 
corresponds to linearizing around u = 0,v = Vq = 2Ae CfJ-x °. Notice that Vq < 0. We get 
then for the coordinate u the following equation (the equation for v is not of concern in this 
context) 

U = CL\V + CL2U . (33) 

It turns out that a\ = 0: this confirms that the trajectory u = is a solution of the system 
(P0|). The sign of determines the stability of the attractor: when a<i is negative, in fact, 
u — > starting from any small perturbation uq. The process evolves with a time-scale 1 / 1 a-2 1 - 
We have 

a 2 = W [(18-/i 2 )/2/i 2 ] (34) 
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which is indeed negative because vq is negative and m < a/3. The stability time scale 1 / 1 «2 1 
goes like (/3)~ 1//2 , i.e. like £A -1 / 2 : the process of infalling toward the attractor is faster as 
£ — > or A — > oo. 
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FIGURE CAPTIONS 



1. The phase-space portrait of the model with f{<p) = </> 2m and V((ft>) <~ </> 2 ™, where n = l,m = 2 and 
negative £ is displayed. Notice the two outward directed attractors and the central Fricdmann region. 
This model is outside the successful range of parameters. All the trajectories spring out from the 
singular points at infinity on the Poincare circle at angles 7r/2 and 37r/2. The diagram is symmetric 
with respect to the origin. 

2. As in Fig. 1, with n = 5, m = 2. This model is inside the successful range: the power-law attractors 
have a natural end. 

3. As in Fig. 1, with n — 4, to = 2. This model satisfies the condition M = 2. Along the deSitter-like 
attractors one has = (f) — ► 0. 

4. As in Fig. 1, with n = 0, m = 1. This phase-space is equivalent to extended inflation. 
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